Semiclassical Theory of Coherence and Decoherence 
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A general semiclassical approach to quantum systems with system-bath interactions is developed. 
We study system decoherence in detail using a coherent state semiclassical wavepacket method 
which avoids singularity issues arising in the usual Green's function approach. We discuss the 
general conditions under which it is approximately correct to discuss quantum decoherence in terms 
of a "dephasing" picture and we derive semiclassical expressions for the phase and phase distribution. 
Remarkably, an effective system wavefunction emerges whose norm measures the decoherence and 
is equivalent to a density matrix formulation. 
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I. INTRODUCTION 

The challenge of understanding to what extent a quan- 
tum system can retain its coherence in the presence 
of interactions with other degrees of freedom has at- 
tracted much attention. Much of it is motivated by 
recent advances in mesoscopicpl 0, Q and cold-atom 
experiments 0, Q as well as keen interest in quan- 
tum computing, which depends crucially on quantum 
coherence 13- 

An essential ingredient in any discussion of coherence 
and decoherence is the identification of a "system" and 
a "bath" which interact with each other in such a way 
that a meaningful distinction can be made between the 
two. In the double slit experiment, for example, the elec- 
trons (or other quantum objeets, e.g., "bucky balls" Q) 
are taken as the system and the degrees of freedom in 
the slits (phonons or spins, for example) are taken as the 
bath. Because the experiment only involves detecting 
the interference pattern on a screen behind the slits, no 
direct measurement of the bath (i.e. the slit degrees of 
freedom) is made. (Since one has no knowledge of the 
state of the bath one must sum over all possible states of 
the bath, i.e. one "traces over the bath".) Only the sys- 
tem is directly observed. As is well knownQ, if the bath 
deteets the path of the particle no interference pattern 
will be seen; the particle has therefore decohered. On the 
other hand, if there is no or only partial detection of the 
path of the particle by the bath, some interference pat- 
tern will be seen with its intensity reflecting the degree 
of coherence of the particle 0- We will later show how 
these familiar statements appear in a very transparent 
way in our semiclassical formalism. 

The problem of a quantum system interacting with an 
environment has been addressed many times in the liter- 
ature, e.g. [lCt llll Il2fl. The Feynman-Vernon influence 
functional approach is well known, although its useful- 
ness beyond the context of harmònic baths has been an 
issue. The influence functional approach to more real- 
istic systems has been advanced significantly by Makri 
and Thompson|l3l ITij . exploiting and developing coher- 



ent state methods with smooth kernels suitable for Montc 
Cario sampling. However, the generality of their ap- 
proach necessarily means some detailed insights and lim- 
iting cases are lost in the machinery, so to speak. 

Another approach is to make a semiclassical approx- 
imation for the system-bath evolution. Casting the 
system-bath interaction and the decoherence problem in 
terms of semiclassical wavepacket dynamics proves to be 
a useful and insightful exercise. Since a semiclassical ap- 
proach is based on classical trajectòries, one can present 
an intuitive picture of the criteria for coherence and de- 
coherence in a coupled system-bath. One may imagine a 
more traditional van Vleck semiclassical Green's function 
approach; however, this has the dimculty that càustic in- 
finities abound in the van Vleck prefactors due ultimately 
to the failure of the stationary phase approximation in 
the limit of small action changes. For example, suppose 
we consider a harmònic oscillator and represent the n th 
state, semiclassically so that (x\n) has singularities 
at the classical turning points for energy E n . Suppose 
now we displacc \n) slightly in position; call this \n). 
The semiclassical projections onto all of the original ba- 
sis states (n\n) are completely wrong for small displace- 
ments. All but one of the projections are incorrectly 
predicted to be since their classical manifolds do not 
overlap, whereas the overlap with the undisplaced orig- 
inal state is nearly singular. The same displacement of 
the n th harmònic oscillator state, expanded in terms of 
localized Gaussians, is quite aceurate; it has no such dif- 
ficultics. Slight displacements of system or bath states 
is commonplace in the decoherence problem, so we avoid 
the càustic difficulties by starting with a wavepacket de- 
scription, avoiding the singularities. We have called this 
the "oil on troubled wàters" effect of using a wavepacket 
description|l5j. 

One of the insights which emerges from this approach 
relates to recent discussions in the literature concerning 
the equivalence of a "bath overlap" picture of decoher- 
ence with a "system dephasing" picture. (See Stern, 
Aharonov and Imry (SAI)|2|, [13, Lüi and Feynman and 
Vernon^çj-) We find that there are three main processes 
that contribute to decoherence: (i) Phase jitter (ii) Bath 
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overlap decay and (iii) Shifts in the trajectory of the sys- 
tem wavepacket. We present explícit formulas for each 
of these effects within our semiclassical wavepacket de- 
scription. 

This paper is organized as follows. In Sec. [H] we de- 
velop the main ideas of our paper and derive the general 
expression for the coherence of a quantum system cou- 
pled to a general bath (described by a density matrix), 
Eq. (|13[1 . This expression can be recast in the form of 
a very intuitive effective system wavefunction, Eq. (|14|l . 
which makes transparent the effects of system-bath in- 
teractions on the system (described by the norm of the 
effective wavefunction). In Sec. IIIII we study several im- 
portant cases of Eq. and Eq. 114(1 in which certain 
physical process (phase jitter, etc.) dominate the sys- 
tem decoherence. In Sec. IIVI we summarize our main 
results and conclusions. Important supplementary mate- 
rial is presented in the appendixes. In Appendix [S] we 
detail how to compute the equations of motion pertur- 
batively for Guassian wavepacket dynamics and derive 
expressions needed in the main text. In Appendix iBl wc 
sketch the arguments of Stern, Aharonov and rmrv|12j 
which equate "bath overlap" and "dephasing" in a spe- 
cial case system-bath interaction. 



II. SEMICLASSICAL THEORY OF 
DECOHERENCE 

We set out to construct a general formal context for de- 
coherence, with the goal of reaching a useful and intuitive 
physical picture. The most general formal structure for 
decoherence (e.g. influcncc functionals for general anhar- 
monic baths) would not involve semiclassical approxima- 
tions, and could claim formal exactness. However such 
formulations must necessarily miss the mark on the issue 
of "useful and intuitive" . 

It is helpful to have a specific model in mind. The 
model of a two armed device already introduced and 
used, for example, in the work of Stern, Ahronov and 
Imry^3 serves that purpose well. In Fig.^ a wavepacket 
representing the system is coherently split into two 
pieces, one of which later interacts with a bath. The 
degree of coherence can be checked experimentally by 
combining the packets (as in the Aharonov-Bohm exper- 
iments of Ref . ]17\ ) . although it is more general to check 
M. = Tr[p^ cd ], where /5 re d is the reduced density matrix 
for the system, after tracing over the bath variables. 

An ambiguity in checking for interference fringes of the 
recombined beams is illustrated by supposing that there 
is a static potential maximum in the left arm but not the 
right, and no system-bath interaction in either arm. This 
will cause a time delay of the left wavepacket compared 
to the right, so no interference will result, even though 
the system is completely coherent. There are ways to 
avoid this, such as taking longer coherent wave trains 
initially for the system, but one must be careful, and 
other ambiguities can arise. 




FIG. 1: A wavepacket representing the "system" in a narrow 
waveguide is split coherently between two arnis, one of which 
interacts with a "bath" of partides. Upon recombining on the 
other side, decoherence due to the system-bath interaction 
affects the interference (or lack of it) between the two arms. 



On the other hand, M. is simply 1 for a completely 
coherent system, and less than 1 if some decoherence has 
taken place. For example, suppose we have a system 
wavepacket broken into two non-overlapping but coher- 
ent pieces, i.e. ip = l/v / 2(V' í +i> r ): with = and 
(ip l \ip l ) = (ip r \i/j r ) = 1. Then the density matrix, p, is 

with i and j taking on the vàlues / and r. It is easily seen 
that Tr[p] = Tr[p 2 ] = 1, i.e. the system is completely 
coherent. However, if somehow the two parts / and r 
become completely decohered, (this in fact requires the 
action of more degrees of freedom-a bath) we loose the 
off diagonal elements of p, getting 

~p = \{W){^ l \ + mm) ■ (2) 

Now we have Tr[p\ = 1, but Tr[,5 2 ] = 1/2. 

This is not the end of decoherence for this system, if 
the left and right packets somehow undergo their own, 
"internal" decoherence. This "internal" decoherence will 
happen on a much longer time scale than the decoher- 
ence of the two initially coherent wavepackets because it 
is the difference in the interactions that each wavepacket 
experiences that determines the decoherence rate. This 
rate is almost always larger for two separated wavepack- 
ets than for a given wavepacket. All this is elementary, 
but it sets the stage for the more detailed work to follow. 

Turning now to a discussion of our general problem 
of a system interacting with a bath, we cast our expres- 
sions in terms of density matrices. The most general 
density matrix for a bath expressed in terms of Guassian 
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wavepackets is 



Pbath 



',w(i,j)\G i0 )(G j0 \ 



(3) 



where w(i,j) satisfies all the properties necessary so that 
/Sbath is a well-defined density matrix (Ti'bath [pbath] = 
1, [u;(i, j)]^ — w*(i,j) — w(j,i), etc). In the special case 
of a pure state bath, we can write w{i,j) = WiWjé 1 ^ - ^ 
where the Wi are real, positive numbers and the fa are the 
associated phases. The states |Gjo) are rnultidimcnsional 
Gaussian wavepackets representing the bath states. Ex- 
panding the bath states in Gaussians allows us to make 
use of very intuitive notions of classical mechanics (which 
guide the wavepacket trajectòries in the semiclassical ap- 
proximation) while at the same time permitting us to 
overcome technical difficulties (singularities) in a more 
naive semiclassical treatment with van Vleck propaga- 
tors. 

Assuming that we take our system be the wavepacket 
initially split into two coherent pieces as shown in Fig. ^ 
the initial wavefunction of the system can be written as 



IWs) = ^= (Wo) + Is5» 



(4) 



where the states \g^) and \g l Q ) are the right (r) and left (1) 
Gaussian wavepackets (coherent states) near the region 
A in Fig. n The density matrix for the system is then 
/5 S ys = |V'sys)(V'sys|- Expressed in terms of Eq. 



the piece of the wavefunction that interacted with the 
bath on the right arm. This approximation actually bc- 
comes exact as h — > 0. It is precisely the lack of "lo- 
cal entanglement" (Eq. Q) still implies "global entan- 
glement" of course) in the wavefunction that makes our 
approach conceptually convenient. In a more general ba- 
sis, we would have a sum of terms for the right arm piece 
of the wavefunction and it would be difficult to identify 
physical phases of the type given by Eq. (|A19(1 . 

The total density matrix at later times thus becomes 

PtotW = ^5>(i,j) {\à t )\G\ t ) + \g r tt )\G r u )e^) 

i, 3 

x((4KG^| + < ff J t |<G; t |e-^08) 

which can be rewritten as 

Ptot(í) = P ll (t) + p lr (t) + f?\t) + jT(t) , (9) 

where p u (í) , p lr (t) , etc. have the obvious meaning. To 
study the coherence of the system, we trace over the bath 
degrees of freedom to obtain the reduced density matrix, 

Picd = Trbath [Ptot] 

= Tr bath [p rr + p H + p w + p u ] , 



which yields, for example, 



(10) 



(ii) 



Av = 2 (Wo) + \9o)) «.9ÓI + (.951) ■ (5) 

The total initial density matrix of the system and bath 
is then 



Ptot (0) — p sys pbath 



(0) 



To study the decoherence of the system due to interac- 
tions with the bath, we must compute the time evolution 
of Eq. Our approach is to use the the perturbative 
wavepacket time evolution described in Appendix^] The 
key result is that the state 



(tò + k r »|G 



,0/ 



g l ot )\G[ t ) + \gl t )\G\ t )e 



Í(pi 



(7) 



where is the wavepacket of the system particle 

moving in the right (left) arm at time t if the bath was 

in state \Gíq) initially. The state \G r i ^) is the perturbed 
wavepacket of the bath (according to Appendix a t 
time t for a particle moving in the right (left) arm where 
the bath was initially in state \Gío). Note that since the 
particle in the left arm does not interact with the bath, 
\ G \t) = \ G it) and \Út) = Wat)- Th e phase & is given by 
Eq. i|A"Ï9)) . 

We emphasize that it is the local nature of the Gaus- 
sian wavepackets combined with weak interactions that 
allows us to write down Eq. Q with a product state for 



where 



-)0r 

3 i 



(Gj t \Gl t ) , 



(12) 



with a similar expression for p l f cd (just the hermitian con- 
jugate of p^ d ) and the other terms. The bath wavepacket 
|G° t ) has a superseript to indicate that it is unper- 
turbed from its trajectory if the system travels in the left 
arm. Note that \g l ot ) does not interact with the bath and 
therefore does not develop an index i depending on the 
bath state. 

When p rec i is squared and traced over to obtain the 
decoherence measure M = Tr [p^ Qd ] , the terms A^coh 
'IV [pre d p l T r ed + P^cdPvíd] contain all the information on 
inter-arm coherence. These give 

■Mcoh = lY,(9h\93tM^J>U^)0% n 0^e^-^' 
iji'j' 
= (\]/ s y s |\i/ s y s ) j 



with 



\* aya ) = ^Y,™(3,i)G%e i <i>i\gr jt 



(13) 



(14) 



Remar kably, A^coh is the self-overlap of a (generally non- 
normalized) effective system wavefunction. The emer- 
gence of a wavefunction form is unexpected because we 
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have not specified that the bath was initially in a pure 
state; it may be in a mixed state such as a thermal bath. 
We can check Eq. 114(1 in the limit of no interaction with 
the bath: then, the overlap factors are all unity, the 
phases 4>j — 0, and all the gaussians \g T j t ) are the same 
(normalized) unperturbed state \gQ t ). Then, 



-^=^w{i,i)\g r ot ) = 



V2 



9ot) 



(15) 



implying (\I/ s y s |iI/ s y s ) = M co h — 1/2, i.e. maximum co- 
herence. 

Eqs. 1(13(1 and 1(14(1 are the central results of this section 
and they are the main formulas of this paper. They apply 
to any bath (harmònic or not, pure state or mixed state) 
with weak system-bath coupling and any number of total 
degrees of freedom. 

Eq. 1(13(1 states that the system coherence is determined 
by an effective system wavefunction with wavepacket 
overlap factors (gj'jffjt), bath overlap factors O J , 
weighting factors (coming from the initial bath condi- 
tions) w(j,i), and finally phase factors e l ^ j . The phases 
4>j are classical actions divided by H and are given by 
Eq. I(A19(I . The "system wavefunction" overlap form, 
Eq. 1(14(1 . is especially convenient for intuition and compu- 
tation. Decoherence shows up as a reduction in the norm 
of the system wavefunction. This comes about from any 
or all of three factors: bath overlap decay, phase jitter, 
and system wavepacket displacement. 

The are several more interesting facets of Eq. 1(13(1 de- 
serving discussion. We will do this systematically, by 
considering important special cases which highlight as- 
pects of this formula. 



Suppose that the system wavepackets |gj t ) are not dis- 
placed by the interaction with the bath, then we have 

<9h\9Ït)™l (16) 

and Eq. ((13(1 becomes 



iji'j' 



- 2 w 



(17) 



If our bath had initially been in a pure state, uu(j, i) 



= / dnx*(v)Xr(r)) 



(18) 



which is just a simple bath overlap. Here Xrü) is the 
state of the bath if the particle went around the right 
(left) arm; rj represents the set of bath coordinates. This 
is the notation of AppendixlBl 

In many physical situations, it may be the case that 
the largest contribution to the sum in Eq. 1(18(1 comes 
from i = j. (This might be the case because the bath 
overlaps are small for i ^ j and/or because the off- 
diagonal terms in w*Wj oscillate in sign from term to 
term.) Making this approximation, we find 



fi 



(19) 



III. SPECIAL CASES OF DECOHERENCE 

We begin our discussion of various limits by assum- 
ing that the system overlaps and possibly also the bath 
overlaps are near unity. This regime is indeed accessible, 
since the classical action perturbation term is cf>i = j-SSt, 
where 5S t is the action due to the perturbation along the 
unperturbed orbit. There is no doubt the action term 
can be large compared to 2tt and vary widely, since the 
perturbing classical action can be large compared to h. 
At the same time, the wavepacket displacement can re- 
main small compared to its width, in both position and 
momentum space. Classical action changes are always 
accompanied by corresponding areas or volumes in phase 
space; if one plots the manifolds of the perturbed system 
exactly, then a phase <f> = 2tt will be accompanied by a 
loop or area in phase space which is of this magnitude. 
However, the wavepacket width goes as ~ y/h, but the 
perturbing action increases as fi. -1 . Thcrefore, for small 
enough perturbations and small enough h, we can safely 
take the wavepacket overlaps to be 1, and focus on the 
phase terms. 



This gives the interpretation of M C oh as one half the 
modulus squared of a phase factor times a bath overlap 
averaged over different "runs" , or realizations of the bath, 
i.e. 



1, 



M coh = -| <e***Crçr > Wi 



(20) 



From the system wavefunction viewpoint, we have (in the 
limit of Eq. IÓ), 



9at) 



(21) 



which lays the blame for decoherence entirely in the sum 
contained in the parentheses. This can be reduced in 
magnitude by both bath overlap decay factors or by 
phase jitter. 



A. Nondynamical Bath 

An important special case to consider is one in which 
the bath does not have any dynamics of its own, the so 
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called "nondynamical" bath. The "nondynamical" bath 
limit emerges by further setting 0®J = 5^ in Eq. I|17|) (or 
0?J = 1 in Eq. lUSt ), i.e. bath wavepackets undisplaced 
by the interaction, which would be the case indeed if 
the bath Hamiltonian commutes with the bath-particle 
interaction potential. Thcn, 




The reduction of the norm (corresponding to decoher- 
ence) is due entirely to phase jitter. In this case we have 
a compelling formula emerging, in the spirit of SAI (See 
Appendix [BJ , 

/ dvxtWxAv) = E W 2e#I ^< e< * > . ( 23 ) 

i 

where the phase 4>i is imparted with probability \uüí\ 2 : 
<pi is the phase acquired if the bath wavepacket is \Gíq) 
initially, and this happens with probability \wí\ 2 , the 
probability weight of that wavepacket in the initial bath. 
This formula gives a concrete picture of the nondynam- 
ical bath limit, and the origin of the phases which are 
averaged over: they are classical actions for the trajec- 
tory of the system-bath dynamics, divided by h. In terms 
of .Mcoh, we have 

Mcoh = Í|<e^>| 2 . (24) 

The limit of a nondynamical bath can also be achieved 
(without the bath Hamiltonian commuting with the 
bath-particle interaction potential) by a high tempera- 
ture bath whose wavepacket description involves mostly 
very excited coherent states. Such coherent states are 
robust against self overlap decay, unless large energy ex- 
change is occurring. This kind of nondynamical bath cor- 
responds to classical states of the radiation field in a large 
cavity with high cnough temperature. It is wcll known 
that this situation is described quantum mechanically in 
terms of excited coherent states of the field oscillators. 
Such baths (or similarly, externally applied fields) tend 
not to contribute to decoherence via a bath overlap de- 
cay mechanism for weak coupling, but rather the sort of 
"dephasing" expressed by Eq. (|2*3|> . 

B. Dynamical Bath 

More generally, we see from Eq. (|17fl that the deco- 
herence (still assuming little system displacement) arises 
from two sources: phase decoherence and amplitude de- 
coherence (due to \0%[\ < 1). The latter is caused by 
the bath wavepackets becoming displaced by interaction 
with the system. In this case, it is less compelling to 
associate (é 1 ^) with Yl Iw^O^le 1 "^ , since the overlap 

i 

factors are not naturally written as integrals over phase 



factors, although one could always do this, however ab- 
sent of physical motivation. This situation corresponds 
to SAI's dynamical bath. The present formulation shows 
a pure phase average picture for this case is somewhat 
forced. Bath wavepacket displacement (and in the next 
section, system wavepacket displacement) thus emerges 
as a restraint on a pure "dephasing" picture of decoher- 
ence. 

The dynamical bath limit would be uninteresting if 
decoherence is always dominated by dephasing. But low 
temperature baths are prime suspects for overlap decay 
to dominate dephasing effects. For example, if there is 
just one bath coherent state, e.g. as in a T — bath 
whose true ground state is describe well by a single mul- 
tidimensional coherent state, the Eq. I|13f) becomes 

-Mcoh^lOÏÏI 2 , (25) 

i.e. the decoherence is entirely caused by bath overlap de- 
cay. This is true even if the system wavepacket is strongly 
displaced, since the system wavepacket simply overlaps 
itself in Eq. (|13f) when there is but a single state in the 
sums. It is therefore possible to decohere from the zero 
temperature initial state or a given single coherent state 
of the bath, due to bath overlap decay; however, at zero 
temperature this requires degeneracies of the bath. 

Here, it is especially clear that a phase average pic- 
ture is not natural for a dynamical bath: in this case 
there is only one "quantum trajectory" so to speak, a sin- 
gle product wavepacket which describes the bath-system 
evolution in the right arm. Eq. I|13f) shows that the ef- 
fect in this case is decoherence due to a displaced bath 
wavepacket. 

C. System Overlap Decoherence 

We now relax the artificial (though possible) condition 
that the system overlap terms are essentially 1, i.e. the 
system wavepackets can be significantly displaced by in- 
teracting with the bath. The concepts of a dynamical and 
nondynamical bath still apply; the question at hand now 
is the further role of system overlap decay in forming the 
coherence measure A'ícoh- We assume a nondynamical 
bath to make the analysis simpler; this means that bath 
wavepacket self-overlaps are unity, and any decoherence 
can be blamed either on random phases (the dephasing 
limit) and/or on system overlap decay. We now investi- 
gate the relative importance of these two. 

The relevant effective system wavefunction is then 

I^^TfEKI 2 ^!^)- ( 26 ) 

v i 

The self overlap of this wavefunction is 

M coh = ÍX>| 2 K-! 2 ^^^Mt) • (27) 
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The phase and overlap contributions are manifest. It 
is not possible to give a general rendition of the rela- 
tive importance of phase and overlap contributions to 
this expression; this will depend on the system and bath 
under consideration. However, the diagonal terms al- 
ways survive, even in the limit of strong kicking of the 
system wavepackets. (We remind the reader that even 
though our analysis was perturbative for the wavepacket 
displacements, in the sense of classical perturbation the- 
ory the displacements can be strong in the quantum 
sense. "Strong" is measured by wavepacket overlap de- 
cay, which can be severe even while classical perturba- 
tions are correctly giving the wavepacket displacements. 
See AppendixEl) Restoring the bath overlap factors for 
a moment, for the diagonal terms we get 

■Mcoh« Í£H 4 I°°T < Í£M 4 ■ (28) 

i i 

Since £K| 2 = 1, 

i 

ÇK|*~l/JVT, (29) 

i 

is an inverse participation ratio, where N is the num- 
ber of participating quantum states describing the bath. 
When N is large, Eq. (|2*9"|) predicts that M co h will be 
vastly smaller than 1/2, effectively meaning the system 
is completely decohered in this limit. This is the limit- 
ing form in the strong system-kicking limit and makes 
physical sense: the broader the distribution of quantum 
states in the initial bath (measured by N in Eq. (|29() , the 
more uncertain the "potential" felt by the system (bath 
has a broad distribution of possible states) and hence the 
greater the decoherence. 

There remains a question, however: could the system 
overlap decay strongly as above without strong phase 
randomization, so that a pure dephasing picture would 
miss it? When one considers, for example, the harmònic 
model, the conclusion is soon reached that for finite tem- 
perature it is not easy to strongly displace the system 
wavepackets randomly without strong phase randomiza- 
tion. 



IV. CONCLUSIONS 

In this paper, we have presented a semiclassi- 
cal, wavepacket-based formalism for decoherence. We 
have limited ourselves to the case of a single system 
wavepacket split initially into two mutually coherent 
pieces, one of which interacts with a bath. We derive 
an expression for the measure of coherence in the sys- 
tem, Eq. 113(1 . which determines the coherence in terms 
of wavepacket overlap factors (ffj/ t |ffj t ), bath overlap fac- 
tors 0®j, weighting factors (coming from the initial bath 
conditions) w(j,i), and finally phase factors e l ^ j . The 



phases 4>j are classical actions divided by h and are given 
by Eq. (IOT) . 

One perhaps surprising and potentially very computa- 
tionally and intuitively useful aspect of our formulation is 
the emergence of an effective system wavefunction, which 
measures the decoherence, Eq. (|14f) : Decoherence shows 
up as a reduction in the norm of the system wavefunc- 
tion. Similar ideas have been introduced in the context 
of a stochastic Schrodinger equation |ïsj . 

After the derivation of the general formulas for the co- 
herence of a quantum system interacting with a bath in 
Sec. |ü] we discuss several special limits of the interac- 
tion. In one limit, discussed in Sec. lIII Al neither system 
nor bath wavepackets are significantly displaced by the 
interaction, but a distribution of phases develop which 
decohere the system. This is limit is naturally described 
as "dephasing" and is appropriate to a number of physical 
situations where the system-bath interactions are quite 
weak. In Sec. lIII Bl a situation is discussed where system 
wavepackets are barely perturbed but bath wavepackets 
are significantly displaced. In this limit, one can still 
force a random phase picture, but the identification with 
an average over random phase factors is more of a math- 
ematical equivalence than a physically motivated idea. 
Finally, if the system wavepacket is strongly perturbed 
by the interaction, as in Sec. IIII Cl a new decoherence 
mechanism sets in: system overlap decay. Such system 
disturbance is hardly rare or unlikely. Strong perturba- 
tion of the system can occur with or without significant 
bath displacement. 

Our perturbation treatment has certain similarities to 
linear response theory. For chaotic systems (expected for 
say a liquid or gaseous bath) it could suffer the same crit- 
icism 19] that the actual magnitude of the perturbation 
for which the formalism is vàlid is unreasonably small. 
However, it might benefit from the same saving graces as 
linear response theory; namely, that ensembles of trajec- 
tòries are better behaved than individual trajectòries. 

The distinction we are making between phase random- 
ization versus overlap decay has long been central within 
the context of spectroscopy of systems embedded in a 
bath (see, e.g. Ref. [20]). The concepts of "dephasing", 
"depopulation" , and "pure dephasing" are traditional in 
spectroscopy. Within the context of exponential decay, 
the relation 



i 2 Li\ ± 2 

is legion, where T 2 is the dephasing time, T\ is the de- 
population time, and T 2 * is the pure dephasing time. 
The translation of concepts into the present discussion 
is: "dephasing" — > "decoherence"; "pure dephasing" — > 
"dephasing" . The time T 2 is typically the time constant 
for decay of the initial wavefunction created by absorp- 
tion of a photon, and is measured from the width of of 
an absorption line. (If we had introduced a tunnel cou- 
pling between the two arms of our model device, we could 
also have had a natural population decay-the probability 
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of being in each arm. This is an interesting subject for 
future study). 

The approach we have taken to decoherence is not 
limitcd to the physical circumstances used here. The 
semiclassical wavepacket-perturbation approach should 
be applicable to wide variety of situations and physical 
measurables including electron decoherence in metalspïl 
I22I |24| and studies of the classical-quantum 

correspondence[25j . We hope to pursue some of these 
in the ncar future. 
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APPENDIX A: COHERENT STATES AND 
GAUSSIAN WAVEPACKETS 

We present a brief review of Gaussian wavepacket dy- 
namics for our approach to the decoherence problem. It 
is well known that the problem of the usual kinetic en- 
ergy operator with a time dependent potential at most 
quadratic in the coordinates is exactly solvable, and is es- 
pecially simple in the case of initial wavefunctions which 
are Gaussian; these remain exactly Gaussian wavepack- 
ets under the time evolution. 

We note that our goals extend far beyond such 
quadratic systems; we will see that a semiclassical ap- 
proximation permits the use of quadratic form dynamics 
in more general contexts. We make use of the so-called 
"thawed guassian approximation" [2(| [53, which em- 
ploys the auxiliary variables Z and Pz , whose dynamics 
are given by the equations below. The "thawed Gaussian 
approximation" allows one to approximate the potential 
locally as quadratic thus taking advantage of the exact- 
ness of Gaussian propagation on quadratic potentials. 

In a multidimensional form, a general Gaussian 
wavepacket is given by 

Í>(q.,t) =exp{-[( q- q t ) T • A t ■ (q - q t ) 

+Pt • (q - qt) + s t }} , (Al) 

where At is an N x N matrix for N coordinates describ- 
ing the stability of the center of the Gaussian wavepacket, 
and q t , pt are TV-dimensional vectors describing the po- 
sition and momentum evolution of the center of the 
wavepacket. We have introduced the more conventional 
wavepacket notation qt = q(í), etc. Let the classi- 
cal Hamiltonian, H = T + V, have the usual Cartesian 



kinetic energy operator and a general time dependent 
potential V smooth at least up to quadratic order in 
the coordinates. The parameters of the Gaussian then 
obeyHJIIi 



dt 
d 
dt 



qt = V„H 



Pt 



A t = -P z Z 



d_ /P Z 
dt l z 



-V"(t) 



s t = Lt+ -Tr[Z • Z- 1 ] , 



(A2) 
(A3) 
(A4) 
(A5) 
(A6) 



where L t is the usual classical Lagrangian. Integrating 
over time, we have 



s t = s + S t + l —Tr[lnZ] 



(A7) 



where V" and m _1 are N-dimensional matrices of mixed 
second derivatives of the Hamiltonian with respect to po- 
sition and momentum coordinates, respectively. That is, 



[V] 



d 2 H 

dqidq-j 



(A8) 



and so forth. St is the usual classical action. Eqs. (JA2J- 
IA6(I holds for a general time dependent V. We focus 
on the stability equations, Eq. I|A5J) . which admit the 
solution 



Pz t 
Z t 



where 



M(í)=fe/' K(t,) dt ' 



(A9) 



(A10) 



and T denotes the time ordering operator, needed be- 
cause 



K(í) 



-V"(í) 

Q 1 



(All) 



does not commute with itself (in general) at different 
times. M(í) is the usual classical stability matrix, where 
Mn = dpt/dpo, etc. 

Consider a narrow (in q) Gaussian wavepacket cen- 
tered on the classical position qo and momentum po- 
Assuming a reasonably smooth potential, let us expand 
around qo up to quadratic terms, arguing that the tails 
of the Gaussian are negligible where the Taylor expan- 
sion starts to break down. We use this quadratic form 
to propagate the packet in the next time instant. Thus 
the Gaussian will propagate in the next instant accord- 
ing to Eqs. i|A2IA6|l . If we agree to move the center 
of the Taylor expansion to the moving mean position of 
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FIG. 2: A Gaussian wavepacket in an anharmonic potential 
is approximately propagated by expanding the potential to 
quadratic order locally around the instantaneous center of 
the Gaussian. 



the wavepacket, qt, then Eqs. I|A2IA6(I will hold, since 
the potential is now, by construction, a time depen- 
dent quadratic form. However, the interpretation has 
changed-the position and momentum parameters qt and 
Pt are now just exactly the usual classical trajectòries on 
the exact, anharmonic potential, but the distortion of the 
wavepacket is governed by the local quadratic expansion of 
the potential-thus keeping the wavepacket Gaussian|2y|. 
We illustrate the idea in Fig. [3 

In general this approximation breaks down after somc 
time due to wavepacket spreading, but that time can be 
put off as long as we please as h — *■ 0, since we can 
take a narrower wavepacket, with position and momen- 
tum uncertainties going as Vh. This delays the spread- 
ing by at least a factor ~ l/|log/ï| in time (for chaotic 
systems) 0, El SUI- 

Since any quantum state (aside from spin states) can 
be built out of Gaussians, we have a full semiclassical 
approach, exact as % — > 0. Each Gaussian is propagated 
with its own optimized time-dependent Hamiltonian. 

The phase St of Eq. (|A7(I is the usual action, taken 
along the guiding trajectory, modified by an extra term 
which takes the place of a Maslov phase. This term 
evolves smoothly in time and therefore is another advan- 
tage of a wavepacket approach as compared to the more 
troublesome eigenfunctions of Hermitian operators. 

There are other approaches based on wavepackets 
which share the smoothing property, such as a full semi- 
classical coherent state propagator [3^ |3|| . This dif- 
fers from the above Gaussian wavepacket (GWP) method 
in that the propagator is optimized for both the initial 
Gaussian and a final Gaussian onto which it is projected. 
It is also the "natural" semiclassical approximation if one 
rctains stationary phase as the defining idea while passing 
to a coherent state basis. The coherent state method is 
usually more accurate at the same value of fi compared 



to GWP, but it suffers several complexities. Foremost 
among these is that complex classical trajectòries come 
in, (albeit for real time) with their attendant analytical 
difficulties including Stokes lincs. Therefore, we choose 
the GWP approach as a best compromise between accu- 
racy and complexity. In the semiclassical limit, fi — > 0, 
there is no accuracy compromise. 

Perturbation theory may seem limited in scope at first 
glance, but we soon realize that strong interaction with 
many bath degrees of freedom leads to such ràpid and 
near complete decoherence that an analysis of the strong 
system-bath coupling regime seems a little post mortem. 
We therefore focus our attention on the weak coupling 
limit. 

In the setup in Fig. ifTjl. we have a Hamiltonian given 
by 

È = H p + È h + V = H + XVi , (A12) 

where H p (Hh) is the Hamiltonian of the particle (bath) 
and V is the coupling between them. Suppose we have 
solved the H = H p + problcm, and now we wish to 
include the effects of AVi, the system-bath interaction. 
In accordance with perturbation theory, the first order 
effect is determined by the extra potential felt by the 
wavepacket as it travels over its old trajectory. Assuming 
this perturbation is smooth as a function of coordinates, 
we include only terms linear in the interaction strength, 
A. We can include perturbations to quadratic order as 
well, at the cost of increased complexity, but we defer 
this for future work. A weak interaction with a bath will 
show up in two ways in the wavepacket given in Eq. JHJ: 
(1) Changes to the guiding trajectory, q t and pt and (2) 
changes to the phase St- 

a. Perturbation of the Guiding Trajectory 

Let qo (t) be the solution of the Ho problem for a par- 
ticular trajectory. The first order perturbed solution we 
take to be q(í) = q (í) + A5q(í), p(í) = p (í) + X6p(t). 
By substituting this into Hamilton's equations, we find 
that <5qt , <5pt obey 

d_(S Pt \ _ ( -\"{t)\(5 Vt \ (i(t)\ 
dt \5q t ) - \ m 1 ) \6qt ) + ^ J ' 

(A13) 

where f(t) = — V^'(qo(í)) is the time dependent forcing 
function. The solution is 

(£) =M( t )(^) + M( t )/M( t r'(t ) )*' 

o 

(A14) 

t 

= M(í) J M(íT 1 ( f( *' } ) dt', (A15) 
o 
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since 5po = 5qo = in the present circumstances, with 

M(í)- 1 =Tt e~J" K(í,) dt ' , (A16) 

and T't forces the reverse order of times; i.e. earliest times 
to the left in the series expansion of M(í) _1 . Note that 
{d/dt)M_{t)- 1 = M(í)- 1 K(í). Thus, the perturbation 
of a general trajectory in classical mechanics generates 
a linearly forced oscillator problem. This is not surpris- 
ing, since the small additional potential creates small new 
forces on the particle near its old trajectory. 

b. Perturbation of the Phase 

To see how the phase of the wavepacket changes vmder 
the perturbation, we examine 

Ss t = 5S t 

= S Í(T-V - AVI) dt 1 

= S [ (T — Vo) dt' — S í AVi dt' . (A17) 



We see the change comes in two parts. The first is the 
changc in the "old" action (it involves the old Vo) due 
to the new trajectory. Assuming the new one has not 
wandered too far, we have 



«5 / (T - V )dt' = X Pt Sq t , 



(A18) 



by the stationary principle of the action (the RHS of 
this equation is not zero since the perturbation causes a 
drift in final position, X5qt-) Thus, the wavepacket phase 
change, to first order, is 

, ò~s t 

^ = T 

= \ 5 JiT-Vo-XVi) dt' 

= X P -^Í - U Vl (q (t')) dt' . (A19) 



APPENDIX B: THE DEPHASING ARGUMENTS 
OF STERN, AHARONOV AND IMRY 

Stern, Aharonov and Imry (SAI) have argued^l lï^ 
that decoherence may be thought of as dcphasing, i.e. 
that a quantum particle acquires a broad distribution of 
possible phases so that its phase is "randomized" and 
thereby loses coherence. SAI argue that decoherence due 
to shifting a bath into an orthogonal state and decoher- 
ence due to a particle acquiring a broad distribution of 
possible phases are equivalent. In this appendix we give 
the skeleton of their argument in their original notation 
to provide a counter point for our discussion in the body 
of this paper. 



SAI consider a quantum particle (coordinate x) moving 
around both arms of an Aharonov-Bohm ring threaded 
by magnètic flux with a bath (coordinate rf) that only 
interacts with the particle in the right arm as shown in 
Fig.H Partides moving around the left arm are assumed 
not to interact with the environment. The initial wave- 
function is taken to be 



ip(t = 0) = [l(x)+r(x)]®xo(v) 



(Bl) 



and corresponds to the particle having just entered the 
ring region (near point A in Fig.^l, but not yet interact- 
ing with the bath. Here l(x) (r(x)) is the initial particle 
wavefunction on the left (right) arm, assumed to be a 
wavepacket, and Xo(jl) is the initial state of the bath, as- 
sumed to be localized in the right arm. SAI then take a 
final wavefunction (near point B in Fig. ^| 



ip(r) = l(x, t) ® xi(v) + r (x, t) <g> Xr(v) 



(B2) 



where Xiiv) {Xr{v)) i s the state of the bath had the par- 
ticle gone around the left (right) arm|38[. SAI distin- 
guish between "dynamical" and "nondynamical" envi- 
ronments, which either do or do not involve non-trivial 
dynamics on their own, respectively. From Eq. (|B2J| . SAI 
obtain the result that the interference term is (taking the 
trace over the environment) 



2Rc 



l*{x)r(x) / dr)xt(v)Xr(v) 



(B3) 



which allows one to intèrpret the reduction of the inter- 
ference (loss of coherence) in terms of a reduction in the 
overlap of the bath states for the two patris around the 
ring[39j. SAI then argue that one can make the identifi- 
cation 



(B4) 



(éf) = / dr)xï(r))Xr(v) 



where 



ai<t>\ = 



(Xole^lxo) and where for nondynamical 



environments the phase angle </> is 



dtV[x{t)]/h 



(B5) 



with x(t) the classical trajectory of the particle around 
the right ring arm|40|. For the case of dynamical envi- 
ronments 



(Xo\e^\xo) = (xo\fe^Io dt ' 9 '^\ X0 ) , 



(B6) 



where Vi(t')=é i6bt '/ h Ve~ lÈbt ' ' /h is the potential in the 
interaction representation and T is the time ordering op- 
erator. A nondynamical environment is distinguished 
from a dynamical one in that the interaction picture op- 
erator V}(í') commutes with itself at different times in 
the nondynamical case. For a dynamical environment it 
is not generally possible 0, ^| to write down a simple 
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relationship such as that expressed in Eq. I|B5|I . However, 
for a nondynamical environmcnt, 

(e*) = J dtPfàe** , (B7) 

where P(0) = IxofaW)! 2 ^- 

The central result of the work of SAI is the equivalence 
expressed in Eq. I|B4(I which states that the reduction of 
Eq. (|B3|) can be viewed as either the shifting of the bath 
states for the two paths around the ring (due to the in- 
teractions in the right arm of the ring) or as the particle 
on the interacting arm being subject to an uncertain po- 
tential resulting in an uncertain phase shift and hence a 
reduction in (é 1 ^) . 



The key approximation of SAI is the assumption that 
the interaction in the right arm is weak enough to ne- 
glect changes in the trajectory of the particle, while still 
allowing the bath to change in response to the presence 
of the particle. This approximation climinatcs entangling 
on the upper arm (the overall wavcfunction is, of course, 
still entangled) and results in the appearance of a direct 
product of bath states in the interference term, Eq. (|B3|I . 
This simplihcation is necessary to reach a "pure dephas- 
ing" expression. In the general situation, the trajectory 
of the particle is in fact altered in the interacting arm, 
which lcads to entangling with the bath, a dephasing ef- 
fect not included in the SAI picture. 
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